New Superhard Phases for 3D Ceo— based Fullerites 
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We have explored new possible phases of 3D Ceo-based fullerites using semiempirical potentials 
and ab-initio density functional methods. We have found three closely related structures - two body 
centered orthorhombic and one body centered cubic - having 52, 56 and 60 tetracoordinated atoms 
per molecule. These 3D polymers result in semiconductors with bulk moduli near 300 GPa, and 
shear moduli around 240 GPa, which make them good candidates for new low density superhard 
materials. 
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Superhard materials are of obvious practical and theo- 
retical interest. Among them, diamond is still the undis- 
puted hardest substance, and will take some beatings. Its 
hardness comes from the stiffness of its tetravalent car- 
bon bonds, which are also the basis of the widely used 
diamond-like coatings. However, natural diamond is rare 
and its industrial production expensive. For this rea- 
son, much experimental and theoretical work has been 
devoted to studying novel potential candidates for hard 
materials. For instance, another carbon-based crystal, 
3D sp 2 carbon , was theoretically predicted to be even 
harder than diamond but never synthesized. Very re- 
cently, a high pressure phase of silicon clathrates was 
studied |2| . Its bulk modulus resulted to be 90 GPa, very 
close to that of the normal silicon cubic crystal, leading to 
the speculation that, if constructed with carbon instead 
of silicon, it could be as hard as diamond. A similar 
idea was proposed for /3-C3N4 in relation to /3-Si3N4 ||] 
but because only small quantities of the former material 
are available, it remains to be proven that the substance 
is indeed harder than diamond. A general theoretical 
approach to study hardness in any of its definitions is 
difficult, if not impossible, to perform. Traditionally, it 
has been linked to a large bulk modulus but it has been 
shown recently to correlate better with the shear modu- 
lus of the material diamond is the leading system in 
both cases, with 443 and 535 GPa, respectively. 

The discovery of Ceo in 1985 as a third form of crys- 
talline carbon Q was the kick-off for one of the most stim- 
ulating fields for research in materials science in the last 
decades. For the first years after its discovery and mas- 
sive production, the known solid state phases of pure Ceo 
were nearly ideal molecular crystals ||. However, from 
the very beginning, there were speculations on the pos- 
sibility of obtaining crystals in which the Cgo molecules 
were covalently bonded to each other, by means of high 
pressures and temperatures |7H|. Since then, several pa- 
pers |To|-|l4|] have appeared reporting a variety of exper- 
imental results in this line (see Ref. pa] for an updated 



review of the field) . 

Nuiiez-Regueiro et al. |H| reported the crystal struc- 
tures of partially covalent forms of Cgo- They observed a 
crystal containing linear chains of polymerized fullerene 
molecules and also two structures, one rhombohedral 
and one orthorhombic, of covalently bonded molecular 
planes (2D Ceo)- The structures are formed by bonding 
through 2+2 cycloaddition reactions between neighbor- 
ing molecules with a carbon hybridization change from 
sp 2 to sp 3 ; chains and planes are bound together by Van 
der Waals interactions. More recently, the realization of 
3D Ceo-based fullerite has been reported ]l7[ . At least 
four crystal structures were found to be compatible with 
the X-ray data: rhombohedral, orthorhombic (pseudo 
tetragonal), tetragonal and simple cubic. Of these, the 
first two are formed by stacking covalent rhombohedral 
and tetragonal planes. Blank and coworkers (Tf|] have re- 
ported the obtention of 3D fullerene solids of extremely 
low compressibility, and very recently a superhard struc- 
ture has been also reported in detail by Chernozatonskii 
et al. fl9[) ; the crystal is body centered orthorhombic, 
quasi fee, and each molecule develops two intermolecu- 
lar bonds with each of its 12 neighbors, presenting cor- 
rugated surfaces. From the theoretical side, Okada et 
al. 0] have reported calculations on a possible structure 
for pseudo-tetragonal (body centered orthorhombic) co- 
valently bonded 3D Cgo, which might be obtained by 
applying a uniaxial pressure of about 20 GPa to pseudo- 
tetragonal polymerized layers of Ceo- This phase would 
be metallic with a high density of states at the Fermi 
level, and the Ceo molecules would be highly corrugated. 

In this work we report new 3D Ceo-based fullerite 
phases which have a direct relation with the most proba- 
ble precursors: polymerized linear chains and orthorhom- 
bic planes, while the molecules preserve their convexity. 
We will show that our proposed structures have low com- 
pressibilities and large shear moduli and therefore are 
excellent candidates for hard materials. 
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FIG. 1. Views of the crystal structure of the proposed hard (56-4) fullerene solid. The hexagonal turrets connecting nearest 
neighbours molecules are clearly seen in d). Tricoordinated atoms are shown in black. 



We have proceeded as follows: an initial search of the 
stable or metastable structures was done using Tersoff 
potentials which arc known to provide a good de- 
scription of carbon fullerencs [^2|,^3), including several 
polymeric phases 24 . Once the potentially interest- 



ing structures were found, these were used as starting 
points for more detailed and accurate ab-initio calcula- 
tions. The structure, stability and electronic properties 
of each phase were obtained, allowing the complete re- 
laxation of cell parameters and atomic coordinates, and 
analyzed; all data we will give here (except for the cal- 
culated vibrational frequencies) are those of the final ab- 
initio results. 

The ab-initio calculations were performed using 
a numerical-atomic-orbital density functional (DFT) 
method described in detail elsewhere [ p5[ . It has been 
already applied to large fullerene molecules and nan- 



otubes Pql , and many other systems p7| . The calcu- 
lations are done using the Generalized Gradient Approx- 
imation for the exchange-correlation, as parametrized 
by Perdew, Burke and Ernzerhof p8| . Core electrons 
are replaced by nonlocal, normconserving pseudopoten- 
tials pjj factorized in the Kleinman-Bylander form |}0| , 
whereas valence electrons are described using linear com- 
binations of pseudo- atomic orbitals. In this work we have 
used a split- valence double-^ basis set, supplemented 
with polarization d orbitals (DZP) pl| . The radial cut- 
offs were 4.2 and 5.0 a.u. for the s and p orbitals respec- 
tively. We have carried out tests using a triple-^ basis 
set, finding only minor changes in the lattice constants 
(less than 0.3 %) and elastic moduli. The charge density 
was expanded in plane waves with a large cutoff of more 
than 150 Ry (the exact values vary slightly, depending on 
the volume). With these parameters we have obtained a 
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bulk modulus for diamond of 430 GPa, in good agree- 
ment with known values Q. The relaxation of the cell 
parameters and atomic coordinates was performed using 
conjugated gradient minimizations at fixed pressure. 

Starting from two experimentally observed ID and 2D 
polymerized body centered orthorhombic Cqq structures, 
we obtain two new metastable 3D phases. These have 
52 and 56 tetracoordinated atoms per molecule, respec- 
tively, and we will refer to them as (52-8) and (56-4). The 
(56-4) structure is obtained by applying uniaxial pressure 
along the stacking direction of a 2D polymer, while for 
the (52-8) one the pressure is applied in the plane perpen- 
dicular to the direction of linear chains. In those phases, 
the individual molecules preserve their convexity as well 
as the bonds between Cgo molecules originally present 
in the chains or planes (4 and 8, respectively). In both 
cases, six new bonds are formed with each of the new 
nearest neighbor molecules in the body centered position. 
Both structures are intimately related to one proposed 
by O'Keeffe, based on purely geometric arguments ||: 
it is possible to construct a body centered cubic crystal 
with Th symmetry, having all carbon atoms with coor- 
dination four, connected by bonds of equal length. We 
have also studied this phase, from which, once all inde- 
pendent degrees of freedom are relaxed, a (60-0) struc- 
ture is obtained. From the electronic point of view, all 
the new structures are semiconductors, with bands gaps 
larger than 2 eV. A general view of the (56-4) structure 
is schematically shown in Fig. [j], where also the few re- 
maining sp 2 atoms can be seen. The (52-8) phase differs 
from the (56-4) in that both the (001) and the (010) 
planes look like the (001) plane in Fig. [l|. For the (60-0) 
structure, the (001), (010) and (001) planes, which are 
equivalent by symmetry, are bonded as in Fig. 0-a. 

Table I shows the relevant structural data and phys- 
ical parameters of the new fullerites p2[| . The specific 
volumes lie between those corresponding to graphite (9 
A 3 /atom) and diamond (5.6 A 3 /atom) but due to the co- 
valent intcrmolccular bonds are much smaller than those 
of Van der Waals Ceo and C70 crystals (~ 11.5 A 3 /atom). 
The length of the bonds between tricoordinated atoms is 
about 1.35 A; between tri- and tetracoordinated atoms 
about 1.48 A; the remaining intramolecular bond lengths 
range from 1.48 to 1.70 A. The bonds forming the hexag- 
onal turrets which connect molecules in the body cen- 
tered position are about 1.57 A long, and those connect- 
ing second neighbor molecules 1.62 and 1.66 A. As ex- 
pected, the bond angles are grouped around the values of 
graphite, diamond, and 90 degrees; the latter are associ- 
ated with the turrets and intermolecular second neighbor 
bonds. For (60-0), where only sp 3 bonds are present, we 
obtain bonding distances of 1.475, 1.548 and 1.56 A. Re- 
garding the intermolecular distances between first neigh- 
bor molecules in the (110) plane (i.e., those 
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FIG. 2. Energy vs. volume curves for the (56-4) (solid 
lines), (52-8) (dotted lines) and (60-0) (dashed lines) struc- 
tures. 

corresponding to the maximum number of covalent bonds 
between Cgo molecules), the typical values are around 8.4 
- 8.5 A, which agree very well with those experimentally 
found by Marquez et al. fl7| (8.40 and 8.80 A). 

The number of sp 3 bonds suggested the interesting pos- 
sibility of low compressibility materials, and in fact our 
structures have bulk and shear moduli around 300 GPa 
and 240 GPa, respectively. They should be compared to 
the values calculated with the same method for diamond, 
430 and 560 GPa. According to the empirical relation 
found by Teter between hardness and shear modulus [Q , 
these values would correspond to a hardness of about 
40 Vickers, which would place them among the ten or 
so hardest known systems, at the lower limit of super- 
hard materials. It is also interesting that the hardness 
expected for these phases comes together with a consid- 
erable lightness of the materials, since the densities are 
about 30 % smaller than that of diamond. 

Energy vs. volume curves for the new phases are shown 
in Fig. g. It can be seen that the (56-4) and (52-8) struc- 
tures have almost the same equilibrium energy and are 
more stable than the (60-0) by about 1 eV. A suggested 
(56-4) — > (60-0) transition is found at about 14 GPa, in 
which the (56-4) structure is compressed along one of its 
axes. We have also recalculated the structure of Okada et 
al. |2(J] finding that the (56-4) crystal is lower in energy 
by almost 2 eV, and has a smaller volume. The structures 
obtained in this work (and that of Ref.20) are metastable 
with respect to the corresponding precursor planes and 
chains, and correspond to local minima in configuration 
space. However, as seen in Fig. 2, they have a large sta- 
bility range, and require quite considerable amounts of 
energy to be dissociated. 



TABLE I. Unit cell parameters (A); specific volumes v 
(A 3 /atom) and Bulk Moduli B (GPa) for the new proposed 
structures. 
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FIG. 3. Vibrational densities of states for the Ceo molecule 
and the (52-8), (56-4) and (60-0) structures. 

It is interesting to note that these phases could be rec- 
ognized by vibrational spectroscopy methods. Fig. ^ 
shows the vibrational DOS for the isolated Cgo molecule, 
and the (52-8), (56-4) and (60-0) structures, calculated 
using the potentials described in Ref. pii| j; the bands 
above 1500 cm -1 , which are related to graphite- like 
vibrations in Cgo (the in-plane optical frequencies of 
graphite lie around 1590 cm -1 ), completely disappear in 
the (60-0) crystal; only four isolated optical bands remain 
in the (56-4) compound, and eight in the 52-8 one. 

In conclusion, we have studied new phases of 3D Ceo- 
bascd fullcritcs which seem to have quite interesting fea- 
tures: their bulk moduli are about 70 % of diamond and 
their shear moduli compare with those of the hardest sys- 
tems known; the phases are extremely compact, although 
the crystals would be light compared to most other hard 
materials; the original convexity of the molecules is pre- 
served - corrugation is a costly energy consuming process; 
and the (52-8) and (54-6) structures are directly related 
to polymerized chains and layers, which seem to be most 
likely precursors for 3D polymerized phases, and could 
probably be obtained at extreme, possibly uniaxial, pres- 
sures. The present result, joined to that of the metallic 
phase reported by Okada et al. [^TJ and the latest experi- 
mental data, suggest an exciting scenario for research on 
new fullerene based materials. 
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